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ABSTRACT 

We have completed a new survey for H I emission for a large, well-defined sample of 154 nearby (z <^ 0.1) 
galaxies with type 1 (broad-line) active galactic nuclei (AGNs). We make use of the extensive database of H I 
and optical parameters, presented in a companion paper, to perform a comprehensive appraisal of the cold gas 
content in active galaxies and to seek new strategies to investigate the global properties of the host galaxies and 
' their relationship to their central black holes. After excluding objects with kinematically anomalous line profiles, 

^ \ which occur with high frequency in the sample, we show that the black hole mass obeys a strong, roughly linear 

, relation with the host galaxy's dynamical mass, calculated by combining the H I line width and the optical size 

of the galaxy. Black hole mass follows a looser, though still highly significant, correlation with the maximum 
rotation velocity of the galaxy, as expected from the known scaling between rotation velocity and central velocity 
■ dispersion. Neither of these H I-based correlations is as tight as the more familiar relations between black hole 

mass and bulge luminosity or velocity dispersion, but they offer the advantage of being insensitive to the glare 
of the nucleus and therefore are promising new tools for probing the host galaxies of both nearby and distant 
AGNs. We present evidence for substantial ongoing black hole growth in the most actively accreting AGNs. In 
these nearby systems, black hole growth appears to be delayed with respect to the assembly of the host galaxy but 
otherwise has left no detectable perturbation to its mass-to-light ratio, as judged from the Tully-Fisher relation, 
O ' or its global gas content. The host galaxies of type 1 AGNs, including those luminous enough to qualify as 

quasars, are generally gas-rich systems, possessing a cold interstellar medium reservoir at least as abundant as 
^ i that in inactive galaxies of the same morphological type. This calls into question current implementations of AGN 

' feedback in models of galaxy formation that predict strong cold gas depletion in unobscured AGNs. 
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sume that the BH-host scaling relations can be directly applied 
to AGNs, of arbitrarily high luminosity and redshift. 

Here we propose a new angle to investigate the relationship 
between BH masses and the host galaxies of AGNs. In normal, 
inactive galaxies, it is well known that the stellar velocity dis- 
persion of the bulge tracks the maximum rotation velocity of the 
disk, roughly as Vm = (V2 — y/3) x ct, (Whitmore et al. 1979). 
Although the theoretical underpinnings of this correlation are 
still murky, and the v„, - cr* relation is not as tight as has been 
claimed (Ferrarese 2002; Baes et al. 2003; Pizzella et al. 2005), 
nevertheless an empirical relation between v„, and ct* does ex- 
ist (Courteau et al. 2007; Ho 2007a), which implies that cr* 
can be estimated from ■;;,„. Now, u,„ can be measured straight- 
forwardly through H I observations for galaxies that are suffi- 
ciently gas-rich. In the absence of a resolved rotation curve, 
Vm can be estimated from a single-dish measurement of the H I 
profile, provided that we have some constraint on the inclina- 
tion angle of the disk. With the rotation velocity in hand, we 
can deduce immediately two physical quantities of interest: the 
total galaxy luminosity through the Tully-Fisher relation (Tully 
& Fisher 1977), and, given an estimate of the size of the disk 
(e.g., from optical imaging), the dynamical mass of the sys- 
tem. The advantages of this approach are clear. Since the H I 
is distributed mostly at large radii (e.g., Broeils & Rhee 1997), 
it should be largely "blind" to the AGN core. This allows us 
to circumvent the problems encountered in the optical, where 
attempts to disentangle the bulge from the active nucleus are 
maximally impacted. Applying this principle. Ho (2007b) has 
recently used the integrated profile of the rotational CO line to 
infer certain properties of high-redshift quasar host galaxies. 

Apart from dynamical constraints, the H I observations yield 
another piece of information of significant interest — a first- 
order measurement of the cold gas content. As the raw ma- 
terial responsible for fueling not only the growth of the BH but 
ultimately also the galaxy, one can hardly think of a more fun- 
damental quantity to ascertain. How does the gas content of 
AGN host galaxies compare with that in inactive galaxies? De- 
pending on one's view on the evolutionary path of AGNs and 
the efficacy of gas removal by AGN feedback, active galaxies 
may be more or less gas-rich than inactive galaxies. Does the 
gas content scale with the degree of nuclear activity? Does it 
vary with the nebular properties of the central source? One of 
the major goals of this paper is to try to address some of these 
basic questions, which to date have been largely unanswered. 

2. THE DATABASE 

The analysis in this paper is based on the database presented 
in the companion paper by Ho et al. (2008), to which the reader 
is referred for full details. In brief, a comprehensive sample 
of 154 nearby type 1 (broad-lined) AGNs was assembled, con- 
sisting of new Arecibc0 observations of 101 sources with z <, 
0.1 1, mostly selected from the Fourth Data Release of the Sloan 
Digital Sky Survey (SDSS; Adelman-McCarthy et al. 2006), 
and a supplementary sample of 53 other nearby sources col- 
lected from the literature. In addition to basic H I properties 
(line fluxes, line widths, and radial velocities), we also assem- 
bled optical data for the AGN (emission-line strengths, line 
widths, line centroid, nuclear luminosity) and the host galaxy 
(image, concentration index, size, axial ratio, total magnitude, 
central stellar velocity dispersion). From this material a number 
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of important physical parameters were derived, including BH 
mass, Eddington ratio, morphological type, inclination angle, 
deprojected rotation velocity, H I mass, dynamical mass, esti- 
mated host galaxy luminosity, and certain rudimentary proper- 
ties of the narrow-line region. Distance-dependent quantities 
were calculated assuming Ho = 70 km s"' Mpc"', il„ = 0.3, and 
f7A = 0.7. 

3. ANALYSIS 

3.1. The v,„ — (7* Relation in Active Galaxies 

To set the stage of using u„, as a surrogate dynamical vari- 
able to investigate BH-host scaling relations, we first examine 
the correlation between v,„ and cr* for the active galaxies in 
our sample. Including all objects that have measurements of 
both quantities, the scatter is discouragingly large (Fig. la). 
Closer inspection, however, reveals that many of the extreme 
outliers correspond to objects that have potentially untrustwor- 
thy deprojected rotation velocities because of uncertain inclina- 
tion corrections, as well as a handful of sources whose central 
stellar velocity dispersions were estimated indirectly from the 
width of the [O II] A3727 emission line (following Greene & 
Ho 2005a). Removing these objects improves the correlation. 
As further discussed in §3.5, a sizable fraction of the objects 
in our sample contain non-classical H I velocity profiles, which 
are either single-peaked, highly asymmetric, or both, indica- 
tive of strongly perturbed or dynamically unrelaxed H I distri- 
butions. If we further remove these cases — a cut that, unfor- 
tunately, drastically reduces the sample to only 40 objects — a 
much cleaner correlation between v,„ and ct* emerges (Fig. \b). 
Because of the small number of objects and their limited dy- 
namic range, the present sample is not well suited to define 
the v,„-a^ relation. Nevertheless, the present sample still re- 
veals a statistically significant correlation between w,„ and ct*; 
the Kendall's r correlation coefficient is r = 0.53, which rejects 
the null hypothesis of no correlation with a probability of P = 
98.2%. We overplot on the figure the v„-(j^, relation obtained 
from the 550 "kinematically normal," inactive spiral galaxies 
from the study of Ho (2007a). The ordinary least-square bisec- 
tor fit for the inactive objects, 

logv,„ = (0.80 ±0.029)logcr* + (0.62 ±0.062), (1) 

provides a reasonably good match to the AGN sample. Accord- 
ing to the Kolmogorov-Smirnov test, the probability of reject- 
ing the null hypothesis that the AGN and non-AGN samples are 
drawn from the same parent population is P = 69.2%. We con- 
clude that the two populations are not significantly different. 
The scatter is still substantial, but recall that the zeropoint of 
the Vm - (T* relation depends on morphological type, such that 
early-type systems have a lower value of u,„/(t* than late-type 
spirals (Um/tr* ~ 1.2-1.4 for E and SO, compared with u„,/cr* 
w 1.6—1.8 for spirals of type Sc and later; see Ho 2007a). That 
the zeropoint for inactive spirals seems to roughly match the 
current sample of AGNs suggests that their host galaxies are 
mostly disk galaxies with modest bulge components, consistent 
with the actual estimated BH masses (median Mbh = 1.6 x 10^ 
Mq) and morphological types (median T = 2.5, corresponding 
to Sab) of the objects. 

Since local AGN host galaxies obey the Mbh - relation, 
and we have just shown that they roughly follow the same 
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Fig. 1 . — Correlation between the central stellar velocity dispersion and the maximum rotation velocity for the active galaxies in our sample. Panel (o) plots all 
available data, including those that have uncertain inclination coirections (open circles) and [O II]-based estimates of the central velocity dispersion (open triangles); 
objects with both rehable rotation velocities and stellar velocity dispersions are plotted as filled circles. No selection has been done with regards to H I profile type. 
Panel (b) excludes all objects with single-peaked and/or asymmetric H I profiles, as well as those with uncertain rotation velocities and [O II]-based estimates of 
the central velocity dispersion. Note the significant reduction in scatter. The dashed line represents the fit to the ti,„ — tr, relation for the 550 "kinematically normaf 
inactive spiral galaxies from the sample of Ho (2007a); see text. 



asymmetric H I profiles (panel ajB. From the Kendall's r test, 
r = 0.58 and P = 99.9%. Overplotted on the figure is the rela- 
tion expected from combining equation ( 1 ) with the Mbh - 
relation for local AGNs obtained by Greene & Ho (2006b), 
log(MBH/M0) = 4.O21og(cr,,/2OOkms-i)-l-7.96: 



log(MBH/M0) = 5.11ogw,„-4.4. 



(2) 



As with Figure 1, we ran a Kolmogorov-Smirnov test to see 
whether the present sample of AGNs is drawn from the same 
population of objects used to define the predicted correlation 
between Mbh and u,„; the probability of rejecting the null hy- 
pothesis of no correlation is P = 68.6%. Again, the two popu- 
lations are not significantly different. 

Nevertheless, the Mbh - relation still contains significant 
scatter, even after omitting the kinematically peculiar objects. 
What is responsible for this? We believe that most of the scat- 
ter is intrinsic. The virial BH mass estimates for nearby AGNs, 
based on broad Ha or H/3, currently have an uncertainty of 
~ 0.3-0.5 dex (Greene & Ho 2006b; Peterson 2007), and er- 
rors associated with the inclination correction applied to the H I 
line widths do not seem to be the dominant source of scatter 
(Fig. 2b). Part of the scatter can be attributed to an effect related 
to Hubble type. As mentioned above and discussed at length in 
Ho (2007a), Um/cr* varies systematically with galaxy concen- 

*For the rest of the v,,, —a^ relation as in inactive galaxies, we expect Mbh to be correlated with d,„ . Not surprisingly, this is confirmed in Figure 2, which 
additionally reemphasizes that a tighter relation results after removing objects with single-peaked and/or paper, we have also excluded two objects from the literature 
sample that have suspiciously low values of v„, (Mrk 359, v„, = 38 km s"' ; Mrk 493, Vm = 24.7 km s"' ). Closer inspection of the original H I data shows that Mrk 359 
has a highly unusual profile consisting of a narrow peak and a very broad base (Springob et al. 2005); moreover, the inclination angle given in Hyperleda, (' = 39.5°, 
seems inconsistent with the nearly face-on appearance of its Digital Sky Survey image (see Appendix in Ho 2007a for a discussion on the uncertainties associated 
with inchnation angles listed in Hyperleda), suggesting that we have seriously underestimated the incHnation coiTection. Mrk 493 may suffer from the same problem; 
the optical (SDSS) image of the object appears to be much more face-on than Hypeiieda's value of i = 45°. 



tration or morphological type, and because our sample contains 
a wide range of morphological types (see Ho et al. 2008), we 
expect the zeropoint of the Mbh - v„, relation to be smeared by 
the actual mixture of galaxy types. This is illustrated in Fig- 
ure 2c, wherein the sample is divided into four broad bins in 
Hubble type. The small subsamples make it difficult to see the 
trend in detail, but comparison of the E-i-SO galaxies with the 
later-type spirals clearly shows that the two subgroups are off- 
set from each other, in the sense expected from the variation of 
Wm/cr* with morphological type. 

Taken at face value, another contribution to the scatter seems 
to be connected to variations in AGN properties, as shown in 
Figures 2d and 2e, where we examine possible dependences on 
broad-line region line width and Eddington ratio. Objects with 
broad Ha or H/3 full- width at half maximum (FWHM) < 2000 
km s"', commonly designated "narrow-line" Seyfert 1 (NLSl) 
galaxies, have a slight tendency to be offset toward lower Mbh 
for a given u,„ (Fig. 2d). Since many NLSl galaxies tend to 
have elevated accretion rates (e.g., Collin & Kawaguchi 2004), 
the weak trend with FWHM may be more clearly discerned by 
dividing the sample according to Eddington ratio. Indeed this 
expectation seems to hold, as shown in Figure 2e, where, within 
the limited statistics, there appears to be a monotonic decrease 
of Mbh/v„, with increasing Lboi/i'Edd- The trend with Eddington 
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Fig. 2. — Correlation between BH mass and the maximum rotation velocity deduced from H I observations. The different panels test for possible dependences on 
(a) H I profile type, (h) inclination angle correction, (c ) Hubble type, (d) FWHM of broad Ha or H/3, and (e) Eddington ratio. In panel (a), objects with uncertain 
inclination corrections were removed; in panels (b)-(e), objects with single and/or asymmetric H I profiles were additionally removed. Panel (/) is the same as 
panel (e), but the different Hubble types have been adjusted to a common zeropoint by shifting Vm ■ The solid line represents the predicted correlation for the whole 
sample, equation (2), given the Mbh-ct* relation of Greene & Ho (2006b) and the Um-cr* relation (eq. 1); the dotted and dashed lines show the zeropoint offset 
(in Mbh) for the objects with Lboi/^Edd < 0.1 and iboi/^Edd > 01> respectively. 



ratio is stronger than that with luminosity alone (not shown). 
Could the apparent variation of MBn/vm with FWHM and 
Lboi/i-Edd be a secondary effect related to the dependence of 
MsH^Vm on morphological type? The most active AGNs in 
the nearby Universe live in moderate-mass, relatively late-type 
disk galaxies (e.g., Heckman et al. 2004; Greene & Ho 2007b); 
these are precisely the hosts of high-Lboi/i'Edd systems, such as 
NLSls. We attempt to remove the Hubble type dependence by 
shifting Vm for the different Hubble type subgroups to a com- 
mon reference zeropoint (defined by the entire sample). As ex- 
pected, the scatter goes down, but, interestingly, the trend with 
LboifLEdd persists (Fig. 2f). 

3.2. Correlation between Mbh and Galaxy Dynamical Mass 

Although our H I observations yield no information on the 
spatial distribution of the gas, we can still obtain a rough esti- 
mate of the characteristic dynamical mass of the galaxy by com- 
bining the deprojected rotation velocity from the H I line width 
with the optical diameter. This approach is justified for the fol- 
lowing reason. The size of the H I disk in spiral galaxies over a 
wide range of Hubble types and luminosities tightly scales with 
the size of the optical disk; within 30%^0%, Dn 1/D25 « 1.7 
(Broeils & Rhee 1997; Noordermeer et al. 2005), where D25 



is the optical isophotal diameter at a surface brightness level 
of iiB = 25 mag arcsec"^. Thus, even if the absolute dynami- 
cal mass may be uncertain because of our ignorance of the size 
of the H I-emitting disk, the relative masses should be reason- 
ably accurate. Following Casertano & Shostak (1980), Ho et 
al. (2008) calculated dynamical masses for our AGN sample, 
which, interestingly, show a fairly strong correlation with BH 
mass (Fig. 3). In particular, BH mass correlates more strongly 
with Mdyn, which scales as £>25U^, than with alone. An ordi- 
nary least-squares bisector fit to the entire sample, after making 
the quality cut as above, gives 

log(MBH/MQ) = (1 .29 ± 0. 1 l)log(Mdy„/MQ) - (6.96 ± 1 . 17), 

(3) 

with an rms scatter of 0.61 dex. The slope is formally steeper 
than unity, but this result should be regarded as provisional con- 
sidering that most of the objects span a relatively narrow range 
of Mdyn. Closer inspection reveals that the zeropoint of the 
Mbh - Mdyn correlation also depends on Hubble type (Fig. 3a). 
At a given Mdyn, early-type galaxies have a higher Mbh than 
late-type galaxies. Relative to the best-fitting line for the whole 
sample, at a fixed log Mbh the shift in logMdyn is +0.19, -0.05, 
-0.10, and -0.59 for galaxies of type E/SO, Sa, Sb, and Sc/Sd, 
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Fig. 3. — Correlation between BH mass and dynamical mass of the host galaxy. Objects with uncertain inclination corrections and that have single and/or 
asymmetric H I profiles are excluded. Panel (a) divides the sample into Hubble types, and panel (b) bins the sample by Eddington ratio, after adjusting the different 
Hubble types to a common zeropoint by shifting Mdy,,. The solid Hne shows the ordinary least-squares bisector fit given in equation (3) for the entire sample; the 
dotted and dashed lines show the zeropoint offset (in Mbh) for the objects with Lboi/i-Edd < 0. 1 and Lboi/i-Edd > 0. 1, respectively. 



respectively. If we adjust the zeropoint of these groups by their 

corresponding offsets (in Mjyn), the rms scatter of the resulting 
correlation decreases to 0.55 dex (Fig. lib). As in the case of 
the Mbh - Vm diagram, we find that at a given Mdyn objects with 
higher Lboi/i^Edd tend to have smaller Mbh- 

3.3. The Tully-Fisher Relation for Active Galaxies 

The empirical correlation between Vm and total galaxy lumi- 
nosity, first introduced by TuUy & Fisher (1977), provides yet 
another avenue to assess potential differences between active 
and inactive galaxies. As described in Ho et al. (2008), approx- 
imate optical host galaxy luminosities can be obtained by com- 
paring the integrated photometry with the nuclear spectroscopy. 
Figure 4 plots the host galaxy absolute magnitude (after remov- 
ing the AGN contribution) converted to the B band, versus the 
maximum rotation velocity; overlaid for comparison is the cor- 
responding Tully-Fisher relation for inactive (spiral) galaxies in 
the Ursa Major cluster (Verheijen 2001). Considering first the 
entire sample (Fig. 4fl), Mb.^osi loosely traces Vm, but the scatter 
is again considerable and the apparent correlation is formally 
statistically insignificant (r = -0.23, P = 93.3%). Pruning, as 
before, the objects with unreliable inchnation corrections and 
kinematically peculiar profiles filters out most of the egregious 
outliers, such that the remaining sample falls mostly within the 
locus of inactive spirals (Figs. Ab and 4c). The formal correla- 
tion between M^^host and Vm is stiU low (r = -0.30) and only 
marginally significant (P = 95.0%), probably because of the 



small sample size and limited dynamic range in Vm (90% of 
the sample has w„, = 200±100 km s"^). Apart from the larger 
scatter exhibited by the active sample, an effect that can be at- 
tributed, at least partly, to its broad mixture of Hubble types, 
the approximate nature in which the host luminosities were es- 
timated, and the sensitivity of the blue bandpass to extinction 
and stellar population variations (e.g., De Rijcke et al. 2007), 
there are no other gross differences between the Tully-Fisher 
relation of active and inactive galaxies. Different Hubble types 
define shghtly offset, parallel sequences (Fig. Ah), an effect that 
has been well-documented in the B band. For a given B-band 
luminosity, an early-type spiral has a larger characteristic rota- 
tion amplitude than a late-type spiral (e.g., Rubin et al. 1985; 
De Rijcke et al. 2007). To separate out this effect, we apphed 
small corrections to the zeropoints (in host galaxy luminosity) 
of each Hubble type bin. The adjusted distribution shows no 
obvious segregation by AGN properties, such as Eddington ra- 
tio (Fig. 4c). 

3.4. Hi Content 

We next examine the H I content of our sample, with special 
emphasis on whether AGN hosts differ in any noticeable way 
from inactive galaxies. In absolute terms, our survey objects are 
quite gas-rich in neutral atomic hydrogen. The H I masses for 
the detected sources range fromMn i « 10^ to 4 x 10'° M© (Ho 
et al. 2008). Taking into account the upper limits for the nonde- 
tections, the Kaplan-Meier product-limit estimator (Feigelson 




Fig. 4. — The B-band TuUy-Fisher relation for the active galaxies in our sample. In panel (a), solid symbols denote rehable deprojected rotation velocities, 
while less reliable values are marked with open symbols; crosses indicate lower limits on v,,, that lack estimates of the inclination angle. In panel (b), the sample is 
divided according to Hubble type (open circles = E/SO; triangles = Sa; squares = Sb; filled circles = Sc/Sd), and panel (c) bins the sample by Eddington ratio (solid, 
open, and star symbols denote objects with Lbol/^Edd ^ 0.1, 0.1-0.3, and > 0.3, respectively), after adjusting the different Hubble types to a common zeropoint by 
shifting A/fi hos(- Objects with uncertain deprojected rotation velocities and single-peaked and/or asymmetric H I profiles were removed from panels (b) and (c). The 
solid line represents the B-band Tully-Fisher relation for inactive spiral galaxies, as derived by Verheijen (2001) for galaxies in the Ursa Major cluster; the dashed 
lines mark the region that has twice the rms scatter of the Ursa Major sample. 



& Nelson 1985) yields a mean of Mh i = (7.0 ±0.66) x lO'^M© 
and a median of Mh i = 5.0 x 10^ Mq. This is almost identical 
to the total H I mass of the Milky Way (5.5 x 10^ Mq; Hart- 
mann & Burton 1997), and comparable to the average value for 
Sb spirals (Roberts & Haynes 1994). The above statistics per- 
tain just to the newly surveyed SDSS sources. The literature 
sample is both heterogeneous and possibly biased because the 
database from which the H I data were compiled, HyperledcQ, 
does not report upper limits. Nonetheless, we confirm that com- 
bining both samples does not significantly alter these statistics. 

The H I content of galaxies varies greatly and systemati- 
cally across the Hubble sequence (Haynes & Giovanelli 1984; 
Roberts & Haynes 1994). Thus, in order to conduct a meaning- 
ful comparison of the H I budget of active and inactive galax- 
ies, we must know their morphological types. Among the 154 
sources in our sample that have H I detections or meaningful 
upper limits, 148 have estimates of both their morphological 
type and host galaxy optical luminosity. Figure 5 shows the 
distribution of H I masses normalized to the Z?-band luminos- 
ity of the host galaxy, subdivided into six bins of Hubble types. 
For reference, we computed the distribution of Mh i/ifi.host for 
13,262 inactive galaxies culled from Hyperleda. These repre- 
sent all galaxies in the database, reported to be current up to the 
end of 2003, that have reliable entries of H I flux, morphologi- 
cal type, and B-band magnitude; known AGNs were excluded, 
as discussed in Ho et al. (2008; Appendix). As previously 
mentioned, Hyperleda does not record upper limits for H I non- 
detections, so these distributions should be viewed strictly as 



upper bounds. Since the majority of the galaxies in Hyperleda 
are relatively bright and nearby, the detection rate of H I among 
the spirals should be very high (^90%; Haynes & Giovanelli 
1984), such that the observed distribution can be regarded as 
being quite close to the true distribution. The situation for early- 
type galaxies, however, is quite different, because the detection 
rate is only --15% for eflipticals (Knapp et al. 1985) and --30% 
for SOs (Wardle & Knapp 1986), and so the distributions shown 
in the figure are highly biased. 

Comparison of the active and inactive distributions reveals 
two interesting points. Among mid- to late-type spirals (Sb and 
later), the frequency distribution of Mh i/ifi.host is roughly the 
same for the two populations. However, for the bulk of the sam- 
ple, which comprise bulge-dominated Sa spirals and SOs, active 
galaxies appear to be more gas-rich than their inactive counter- 
parts. The most dramatic manifestation of this effect shows up 
among the ellipticals, although here we are handicapped some- 
what by the small number of sources (9) in the active sample. 
Note that among the E and SO subgroups, the difference be- 
tween the active and inactive samples is far greater than por- 
trayed (although we have no rigorous statistical way to quan- 
tify this) because, as mentioned above, the inactive distribution 
is highly biased by the omission of upper limits. (As noted, the 
literature data for the active sample are likely biased as well, but 
we verified that our conclusions are essentially unaffected if we 
exclude these data from the AGN sample.) The host galaxies of 
nearby type 1 AGNs across all Hubble types are at least as gas- 
rich as inactive galaxies, and among early-type systems their 
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Fig. 5. — The distribution of H I masses normalized to tlie B-band lumi- 
nosity of the host galaxy, as a function of its Hubble type. The black his- 
tograms show the distributions for 13,262 inactive galaxies with H I detections 
listed in Hyperleda. The AGN sample from our study is plotted in red his- 
tograms, with upper limits indicated by dotted lines. The entire sample of 154 
objects is plotted in the penultimate panel from the bottom, while the top pan- 
els show the sample sorted by Hubble type. The number of inactive galaxies in 
each group is labeled, with the number of active objects shown in parentheses. 
The bottom-most panel isolates the quasars. The red histograms correspond 
to objects from the cun'ent H I sample; the blue histograms con'espond to 
PG quasars observed in CO (Ho 2005a); the green histograms correspond to 
z < 0.5 PG quasars with dust masses (Haas et al. 2003); see text for details. 

gas content appears to be markedly enhanced. 



3.5. Environmental Effects 

Crude information on the spatial distribution or dynamical 
state of the gas can be ascertained from the degree to which 
the observed H I line profiles deviate from the classical double- 
horned signature of an optically thin rotating disk. Experience 
with nearby galaxies indicates that tidally disturbed systems of- 
ten exhibit asymmetric, single-peaked, or otherwise highly ir- 
regular line profiles (e.g., Gallagher et al. 1981; Haynes et al. 
1998). Given the modest signal-to-noise ratio of our data. Ho 
et al. (2008) decided not to implement a rigorous scheme to 
quantify the line morphology, but rather adopted a qualitative 
classification scheme in which obviously non-classical profiles 
were simply labeled as "A" (asymmetric), "S" single-peaked, 
or "AS" (combination of both). Among the 66 objects detected 
in our new survey, eight (12%) are classified as "A," five (8%) 
are classified as "S," and 16 (24%) are classified as "AS," for 
an overall frequency of non-classical profiles of 44%. Ho et 
al. (2008) performed a systematic census of physical neighbors 
within a search radius of T.5 around each object. There seems 
to be no clear association between profile pecuUarity and the 



presence of nearby neighbors. While some objects with single- 
peaked and/or asymmetric profiles have plausible nearby com- 
panions, many do not; at the same time, a number of objects 
with apparent companions exhibit seemingly regular line pro- 
files. The presence or absence of kinematic irregularity also 
appears to be uncorrected with any of the global or AGN pa- 
rameters that we have at our disposal. We searched for possible 
differences in the following quantities, but found none that was 
statistically significant: morphology type, galaxy luminosity, 
total and relative H I mass, AGN luminosity, broad-line region 
FWHM, BH mass, and Eddington ratio. 

3.6. Connections with AGN Properties 

The availability of optical data affords us an opportunity to 
explore possible connections between the H I and AGN prop- 
erties of our sample. Concentrating on the SDSS objects, for 
which we have homogeneous optical spectroscopic measure- 
ments, we find that their narrow emission-line ratios place the 
majority of them in the territory of Seyfert nuclei (Fig. 6). 
This is not surprising, since most of the SDSS-selected type 1 
AGNs tend to have relatively high accretion rates (Greene & Ho 
2007b), which generally correspond to high-ionization sources 
(Ho 2004b). The objects detected in H I do not stand out in 
any noticeable way from the nondetections. The same holds 
for the electron densities of the narrow-line region, as can be 
inferred from the line ratio [S II] A6716/[S II] A6731: the H I 
detections are statistically indistinguishable from the H I nonde- 
tections (Fig. 7). Next, we searched for a possible dependence 
of AGN luminosity or Eddington ratio on H I content, either in 
absolute (Mh i) or relative (Mh i/i-B.host) terms, but again found 
none (Fig. 8). There is, at best, a mild trend of increasing Ha 
luminosity with increasing H I mass (Fig. 8fl), but given the 
mutual dependence of the two quantities on distance, we regard 
this result as highly suspect. 

4. DISCUSSION 

4. 1 . Black Hole-Host Galaxy Scaling Relations 

The relative ease with which integrated H I emission can be 
detected in nearby AGNs opens up the possibility of using the 
global H I line width as a new dynamical variable to investi- 
gate the scaling between BH mass and host galaxy potential. A 
quantity that can be readily measured from single-dish spectra, 
the H I line width has been widely used in a variety of extra- 
galactic contexts as an effective shortcut to estimate v„,, the 
maximum rotation velocity of galaxy disks. If v,,, correlates 
with bulge velocity dispersion, a,, as originally suggested by 
Whitmore et al. (1979), can be used in place of cr, to de- 
fine a Mbh - Vm relation to substitute for the more traditional 
Mbh - c* relation. A number of authors have suggested that 
galaxies over a wide range in morphological types obey a tight 
correlation between v,,, and (T», to the point that v„, can actu- 
ally replace, and perhaps should be regarded as more funda- 
mental than, (T, (Ferrarese 2002; Baes et al. 2003; Pizzella et 
al. 2005). This result has been challenged by Courteau et al. 
(2007) and Ho (2007a), who found, using much larger samples, 
that the ratio v„,/a^, shows significant intrinsic scatter and sys- 
tematic variation across the Hubble sequence. Nevertheless, a 
Vm - 0", relation does exist, at least statistically, and in circum- 
stances when (T, is difficult or impossible to measure (e.g., in 
very bright or very distant AGNs), Vm may offer the best or pos- 
sibly the only means of constraining the host galaxy. Such was 
the motivation behind the recent study of Ho (2007b; see also 
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Fig. 6. — The location of the H I sample on the line intensity ratio diagrams (a) 
log [S II] AA6716, 6731/Ha. Objects detected and undetected in H I are plotted as 
open squares AGNs from the Palomar survey of nearby galaxies (Ho et al. 1997a) 
variation with H I detection. 
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solid circles and triangles, respectively. The crosses represent H II nuclei and the 
. The H I sample mostly occupies the locus of Seyfert galaxies, with no apparent 



Shields et al. 2006), who used the width of the rotational CO 
Une, locally calibrated against H I, to infer the masses of the 
host galaxies of high-redshift quasars. The line width method, 
however, can prove to be useful even under less extreme condi- 
tions. As discussed at length in Greene & Ho (2006a), a variety 
of factors conspire to make measurement of cr* very challeng- 
ing in galaxies containing bright AGNs, regardless of their red- 
shift. To bypass this difficulty, many studies resort to using gas 
velocity dispersions measured from narrow nebular lines, but 
this shortcut has its own set of comphcations (Greene & Ho 
2005a). Others skip the kinematical route altogether and in- 
stead use the host galaxy (or bulge, if available) luminosity as 
the variable to relate to the BH mass. However, cleanly 
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Fig. 7.— Distribution of the line intensity ratio [S II] A6716/[S II] A6731, 
which is inversely proportional to the electron density, for objects detected 
(hashed histograms) and imdetected (open, dashed histograms) in H I. 



decomposing the underlying galaxy, much less its bulge, in the 
presence of a bright AGN core is often a nontrivial task, even 
under the most ideal conditions (e.g., Kim et al. 2007). 

Within this backdrop, we explored the correlation between 
BH mass and Vm- As expected, a loose correlation exists be- 
tween these two quantities, but the scatter is enormous, 1 
dex (Fig. 2a). Many of the extreme outhers correspond to ob- 
jects with single-peaked and/or clearly asymmetric H I profiles, 
while others have especially questionable inclination correc- 
tions; removing these produces a cleaner correlation. Note that 
apart from being inefficient (roughly 40% of the sample was re- 
jected), this step imposes no serious complication, since classi- 
cal double-horned line profiles are easy to recognize. Nonethe- 
less, even after this filtering, the scatter is still substantial (0.9 
dex; Fig. 2b). We identified a number of parameters associated 
with the AGN that seem to contribute to the scatter (FWHM 
of the broad emission lines, nuclear luminosity, and Eddington 
ratio), but an important contribution comes from the morpho- 
logical type variation within the sample (Fig. 2c). As discussed 
in Courteau et al. (2007) and Ho (2007a), the zeropoint of the 
Vm—(J* relation varies systematically with galaxy concentration 
(which is loosely related to galaxy morphology or bulge-to-disk 
ratio). If we arbitrarily shift the zeropoints of the different Hub- 
ble type bins to the average relation for all Hubble types, the 
scatter decreases somewhat to ^ 0.7 dex (Fig. 2f). 

More promising still seems to be the correlation between BH 
mass and galaxy dynamical mass, which requires knowledge 
of one additional parameter, namely the galaxy's optical diam- 
eter. Retaining, as before, only the objects with robust incli- 
nation corrections and kinematically normal H I profiles, the 
MBH-A/dyn has an rms scatter of 0.61 dex, which further im- 
proves to 0.55 dex after shifting the different Hubble types to a 
common reference point (Fig. 3). While this scatter is larger 
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Fig. 8. — The variation of nuclear Ha luminosity (panels a and b) and Eddington ratio (panels c and d) with H I content. Objects detected in H I are plotted as 
filled symbols, while the nondetections are marked as open symbols with arrows. 



than that of the local AGN Mbh-ct* relation (0.4 dex; Greene & 
Ho 2006b), our estimates of the dynamical masses leave much 
room for improvement, both in terms of getting higher signal- 
to-noise ratio Une profiles and deeper optical images, which 
will allow more accurate measurements of inclination angles 
and isophotal diameters and better estimates of morphological 
types. 



4.2. The Chicken or the Egg 

One of the most intriguing results from our analysis is the 
tentative detection, both in the Mbh - v„ (Fig. 2/) and Mbh - 
Miyn (Fig. 3i)) diagrams, of differential growth between the cen- 
tral BH and the host galaxy. For a given host galaxy potential 
{Vm or Mdyn), AGNs with higher accretion rates (Eddington ra- 
tios) have systematically less massive BHs, implying that in 
these systems the central BH, still vigorously accreting near its 
maximum rate, has yet to reach its final mass, which, from ob- 
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servations of inactive systems, we know to be a well-defined 
constant fraction of the total bulge luminosity or mass (Kor- 
mendy & Richstone 1995; Magorrian et al. 1998; Haring & 
Rix 2004). Many (but not all) of the objects with still- growing 
BHs turn out to have narrower broad Balmer lines (FWHM <J 

2000 km s~' ) because NLS 1 galaxies tend to have higher accre- 
tion rates. At a fixed Vm, the AGNs in our sample with Z-boi/i-Edd 
> 0. 1 on average have BHs 0. 19 dex (factor of 1 .5) less massive 
than those with Lboi/i-Edd < 0. 1 ; in terms of fixed M^yn, the dif- 
ference in BH masses between the low- and high-accretion rate 
subgroups is ~0.41 dex (factor of 2.6). In their study of type 1 
AGNs with stellar velocity dispersion measurements, Greene & 
Ho (2006a) also noticed that the best-fitting Mbh - relation 
for AGNs has a small zeropoint offset of -0.17 dex relative to 
Tremaine et al.'s (2002) fit of the Mbh -o"* relation for inactive 
galaxies. A result similar to that of Greene & Ho (2006a) has 
been reported by Shen et al. (2008), who additionally note that 
the amplitude of the Mbh - cr* relation for AGNs depends on 
Eddington ratio, in the same sense that we find in our study. 
The overall qualitative agreement among the three independent 
studies lends credence to the idea that the most highly accreting 
BHs are still actively growing. 

Mathur et al. (2001; see also Grupe & Mathur 2004; Mathur 
& Grupe 2005), using the width of the [O HI] A5007 Une as a 
substitute for cr*, a shortcut that allows large numbers of AGNs 
to be placed on the Mbh-o"* relation, proposed that NLS Is 
contain BHs that are undermassive with respect to their bulges, 
by as much as an order of magnitude. Employing the same 
technique, Bian & Zhao (2004) arrived at a similar result, as 
did Wandel (2002), who used bulge luminosities as a reference. 
The claim that NLS Is are "young" AGNs, however, has not 
been universally embraced, in part because different authors, 
using the same methods, have arrived at conflicting conclusions 
(e.g., Wang & Lu 2001), because of worries concerning the re- 
liability of the [O III]-based velocity dispersions (Botte et al. 
2004; Greene & Ho 2005a), and, most seriously, because di- 
rect measurements of stellar a, in NLS Is find no compelUng 
evidence that they deviate from the Mbh - ct* relation (Barth 
et al. 2005; Botte et al. 2005; Greene & Ho 2006b). Part of 
the confusion stems from the very definition of the NLSl class. 
Strictly speaking, NLS Is are classified on the basis of optical 
spectroscopic criteria (e.g., Pogge 2000), which, apart from the 
(fairly arbitrary) line widfli limit of FWHM < 2000 km s"\ 
are neither rigorously defined nor universally followed. That 
many cataloged and well-studied NLS Is have high accretion 
rates (e.g.. Pounds et al. 1995; Collin & Kawaguchi 2004) is 
a statistical reflection of the fact that bright AGNs in the local 
Universe tend to have moderately low-mass BHs (Heckman et 
al. 2004; Greene & Ho 2007b), and that soft X-ray selection 
preferentially favors systems in their high-state (Greene & Ho 
2007a). There is no requirement that NLS Is must have high 
accretion rates. The present sample, in fact, provides two use- 
ful examples. By their line widths alone, both NGC 4051 and 
NGC 4395 qualify as NLS Is, and yet the former has Lboi/^Edd = 
6.3 X 10"^ while the latter, with Lboi/i^Edd = 2.7 x 10"^, formally 
has the lowest Eddington ratio in the sample. 

That Msn/vm and MBn/Mdyn decrease systematically with in- 
creasing Lboi/i-Edd suggests not only that the growth of the BH 
and the growth of the galaxy are not finely synchronized, but 
also that there is a preferential sequence to their coevolution: 
during any particular connmon growth event, star formation pre- 
cedes BH growth. Thus, at least for the local sample under 



consideration, we have a tentative solution to the proverbial 
"the chicken or the egg problem." Phenomenologically, star- 
burst activity comes first, followed by AGN activity — not the 
other way around, and probably not concurrently. This resem- 
bles the scenario envisioned by Ho (2005a, 2005b; see also 
Kim et al. 2006) in his analysis of the star formation rate in 
nearby quasars. Interestingly, the sequence of events seems to 
be reversed at high redshifts, where the best evidence to date 
suggests that galaxy growth actually lags BH growth (Peng et 
al. 2006a, 2006b; Shields et al. 2006; Ho 2007b). One might 
speculate that the pathway for BH growth and its interdepen- 
dence on the host galaxy could be radically different for lumi- 
nous quasars in massive, early-type galaxies, compared to the 
much more modest AGNs under consideration here, which re- 
side mostly in later type, disk systems (e.g., Hopkins & Hem- 
quist 2006). On the other hand, we note that moderate-redshift 
(z = 0.36) AGNs, mostly Seyfert galaxies not too dissimilar 
from those in our sample, also follow the trend of high-redshift 
quasars in showing a higher ratio of BH mass to host galaxy 
(bulge) mass (Woo et al. 2006; Treu et al. 2007). We cannot, at 
the moment, reconcile these conflicting results, which urgently 
need to be verified. 

We end the discussion on BH growth with some cautionary 
remarks. Throughout this study we take at face value that the 
virial BH masses contain no major sources of systematic error, 
other than the ^^0.5 dex uncertainty in the zeropoint (Greene 
& Ho 2005b; Peterson 2007). As discussed in Greene & Ho 
(2007c; Appendix), however, our assumption that the virial 
mass formahsm is invariant with respect to AGN properties is, 
at the moment, more an article of faith than an indisputable fact. 
The size-luminosity relation on which the virial formalism de- 
pends continues to undergo revision and is especially poorly 
constrained at the low-luminosity end (Bentz et al. 2006). The 
very structure of the broad- Une region, and hence the geometric 
("/") factor of the virial relation, may depend on accretion rate 
or Eddington ratio, but presently we have no reliable means to 
quantify this effect (Greene & Ho 2007c). Strictly speaking, 
therefore, the trend of decreasing Msu/vm and MBn/Mdyn with 
increasing Lboi/^'Edd can be interpreted as evidence that Mbh is 
systematically underestimated in high-Lboi/i-Edd AGNs. Resolv- 
ing this degeneracy lies beyond the scope of this work. At the 
moment we also cannot prove that the identified trends are not 
the result of subtle mass-dependent selection effects. In their 
investigation of the local BH mass function of nearby AGNs, 
Greene & Ho (2007b) note that low-mass BHs, which even at 
their Eddington limit are still relatively faint AGNs, preferen- 
tially make it into a magnitude-limited survey such as SDSS if 
they reside in more luminous (more massive) galaxies. 

4.3. Host Galaxies and Environment 

Even though most sizable nearby galaxies harbor central 
BHs, they exhibit very feeble traces of nuclear activity (Ho et 
al. 1997b). Only a tiny fraction of the local galaxy popula- 
tion accrete at rates high enough to qualify them as respectable 
AGNs (Ho 2004b; Greene & Ho 2007b). Why is this? The low 
space density of luminous AGNs in the local Universe surely 
reflects some combination of the overall reduced gas supply in 
present-day galaxies as well as their more placid dynamical en- 
vironments, but precisely which physical parameters — internal 
or external to the system — actual trigger the onset of nuclear ac- 
tivity in any given galaxy remains largely an unsolved problem 
(Ho et al. 1997c, 2003; Ho 2004b). In a similar vein, and espe- 
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cially given the recent discussion of the purported link between 
BH and galaxy growth, it is of interest to know what impact 
AGN activity actually has on the properties of the host galaxy. 

Our new H I survey offers an opportunity to examine these 
issues from some fresh angles. Despite the limited quality 
of the data, our visual examination the optical images do not 
give the impression that the galaxies or their immediate envi- 
ronments are particularly unusual. Some clear cases of inter- 
acting or binary systems exist within a separation of 50 kpc 
(Ho et al. 2008; see Fig. 3g, SDSS J165601.61+211241.2; 
Fig. 4b, SDSS J130241.53+040738.6), but overall the optical 
morphologies of the hosts appear undisturbed and large (i.e. 
nearly equal-mass) physical companions are rare. There are 
also no discernible differences between the objects detected and 
undetected in H I (cf. Fig. 3 vs. Fig. 4 in Ho et al. 2008). 

Excluding the objects with anomalous H I profiles, the rest 
of the sample is not just morphologically normal but also dy- 
namically normal. We draw this inference from our analysis 
of the Tully-Fisher relation for the active galaxies, which, be- 
sides the anticipated larger scatter due to various measurement 
uncertainties, reveals no striking differences compared to reg- 
ular disk galaxies (Fig. 4). As an aside, it is worth making an 
obvious point: the fact that our objects obey the Tully-Fisher re- 
lation implies that the H I must be roughly regularly distributed, 
at least regular enough so that its integrated line width traces the 
flat part of the galaxy rotation curve. Hutchings et al. (1987) 
and Hutchings (1989) previously failed to see a clear Tully- 
Fisher relation for their sample of AGNs because of the high 
frequency of asymmetric H I profiles. By contrast, the more ex- 
tensive study by Whittle (1992) found that Seyfert galaxies do 
define a Tully-Fisher relation, but one that is offset toward lower 
velocities compared to normal spiral galaxies. He interpreted 
this as an indication that Seyfert galaxies have lower mass-to- 
light ratios (by a factor of 1.5- 2), possibly as a result of en- 
hanced star formation. Our results do not agree with Whittle's. 
We have verified that the discrepancy between the two studies 
does not stem from the different fiducial relations chosen for in- 
active spirals. While Whittle's reference Tully-Fisher relation 
has a much steeper slope than the one we adopted (from Ver- 
heijen 2001), the two zeropoints are very similar. We speculate, 
but cannot prove, that the low-velocity offset seen by Whittle 
may be caused by contamination from objects with kinemat- 
icaUy peculiar, preferentially narrow H I profiles. Inspection 
of our data (Fig. 4) reveals that prior to excluding the objects 
with asymmetric or single-peaked profiles our sample also has 
a mild excess of low-velocity points. Indeed, a nonnegligible 
fraction of all galaxies, irrespective of their level of nuclear ac- 
tivity, show this behavior, which Ho (2007a) attributes to dy- 
namically unrelaxed gas acquired either through minor mergers 
or primordial accretion. 

Enhanced star formation frequently seems to precede nuclear 
activity (Kauffmann et al. 2003; see discussion in Ho 2005b), 
but at least for the AGN luminosities probed in this study, the 
extent of the young stellar population has left no visible im- 
print on the mass-to-light ratio of host galaxies, insofar as we 
can gauge from the Tully-Fisher relation. 

Do AGN host galaxies possess a higher fraction of kinemati- 
cally anomalous H I profiles than inactive galaxies? Our study, 
as in previous ones (Mirabel & Wilson 1984; Hutchings et al. 
1987), tentatively suggests that the answer is yes, but a quan- 
titative comparison is difficult. Single-dish surveys of isolated, 
(mostly) inactive disk galaxies also report frequent detections 



of asymmetric H I line profiles (e.g., Baldwin et al. 1980; Lewis 
et al. 1985; Matthews et al. 1998), which are attributed to non- 
circular motions, unresolved companions, or genuine perturba- 
tions in the spatial distribution of the gas. Richter & Sancisi 
(1994) and Haynes et al. (1998), for example, find asymmet- 
ric profiles in roughly 50% of the spiral galaxies they surveyed. 
Similarly, Lewis (1987) studied a large sample of face-on spiral 
galaxies and found narrow H I profiles to be rare. He interprets 
this result to imply that a substantial fraction of the galaxies 
must have distorted H I distributions, whose large-scale motion 
is misaligned with respect to the plane of the stellar disk. Thus, 
in absolute terms the frequency of non-classical H I profiles in 
our sample (44%; §3.5) is very similar to the frequency reported 
for inactive galaxies, but we are wary to draw any firm conclu- 
sions from this because the signal-to-noise ratio of our spectra 
is generally much lower than those of the control samples, and 
because of the qualitative nature of our profile classification. 

The kinematically anomalous objects account for 17% of 
the 792 galaxies studied by Ho (2007a), but this fraction de- 
pends on Hubble type, being more common among early-type 
systems. Among SO galaxies, the fraction reaches 40% (Ho 
2007a), nearly identical to the frequency found among the ac- 
tive systems studied here, which, although strictly not all SO 
galaxies, nonetheless tend to be bulge-dominated disk galaxies. 
As in the case of inactive galaxies (Ho 2007a), the kinematic 
peculiarity for the AGN sample cannot be attributed to tidal in- 
teractions or comparable- sized nearby neighbors. 

4.4. Gas Content and Implications for AGN Feedback Models 

Insofar as their H I content is concerned, the host galaxies 
of the AGNs in our sample are endowed with plenty of gas. 
The most meaningful metric is the specific rather than the ab- 
solute gas mass, and because this quantity spans a wide range 
across the Hubble sequence (Roberts & Haynes 1994), we can 
sharpen the comparison even further by specifying the morpho- 
logical type of the galaxy. This exercise (Fig. 5) convinces us 
that type 1 AGNs are at least as gas-rich as inactive galaxies, 
and among early-type hosts, beginning with Sa spirals and cer- 
tainly by the time we reach SOs and Es, AGN hosts appear to 
be even more gas-rich than their inactive counterparts. For ex- 
ample, the sensitive survey of 12 inactive E and SO galaxies 
by Morganti et al. (2006) detected H I emission in nine ob- 
jects, of which only one has Mh i/Z^B.host > 0.1. By contrast, 
among the 30 active E and SO galaxies in our survey that were 
detected in H I, 26 (87%) have Mh i/Lfi.host > 0.1. Our results 
are quahtatively consistent with those of Bieging & Biermann 
(1983) and Mirabel & Wilson (1984). Because of the nature 
of the control sample of inactive objects (only detections are 
available), we hesitate to be more specific, but E and SO galax- 
ies hosting AGNs may be overabundant in H I by as much as 
an order of magnitude. Previous studies of nearby elliptical 
and SO galaxies have hinted of a possible statistical connection 
between H I-richness and radio activity (Dressel et al. 1982; 
Jenkins 1983; Morganti et al. 2006), two well-known exam- 
ples being the low-power radio galaxies NGC 1052 (Knapp et 
al. 1978) and NGC 4278 (Raimond et al. 1981). The validity 
of our conclusion, of course, critically depends on the accu- 
racy of the morphological types. However, to bring the aver- 
age Mh i/i-B.host ratios of the E and SO subgroups to conform 
to the average values of inactive galaxies, we would have to 
shift the morphological types to as late as Sb or even Sc spirals. 
It seems unlikely that the morphologies could be so blatantly 
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wrong. Alternatively, perhaps we have grossly miscalculated 
the host galaxy optical luminosities. Without improved optical 
imaging (§5), however, these alternative explanations remain 
purely hypothetical. 

The global gas content bears no relationship to the level of 
AGN activity, either in terms of absolute luminosity or Edding- 
ton ratio (Fig. 8). Ho et al. (2003) had come to the same conclu- 
sion for nearby LINERs and Seyferts, but now the same holds 
for AGNs on average 2-3 orders of magnitude more luminous. 
Whether H I is detected or not also makes no impact whatso- 
ever on the optical spectrum (Figs. 6 and 7). Perhaps none of 
this should come as a surprise. The H I data, after all, probe 
spatial scales vastly disproportionate to those relevant for the 
AGN central engine. Moreover, the accretion rates required to 
power the observed activity are quite modest. The Ha lumi- 
nosities of our sample range from ~ 10"*'^ to 10** ergs s"\ with 
a median value of 4 x 10^^ ergs s~'. Adopting the Ha bolo- 
metric correction given in Greene & Ho (2007b) and a radiative 
efficiency of 0.1, this corresponds to a mass accretion rate of 
merely ~ 2 x 10"^ M© yr~'. 

Many current models of galaxy formation envision AGN 
feedback to play a pivotal role in controlling the joint evolution 
of central BHs and their host galaxies. During the major merger 
of two gas-rich galaxies, each initially seeded with its own BH, 
gravitational torques drive a large fraction of the cold gas to- 
ward the central region of the resulting merger remnant. Most 
of the gas forms stars with high efficiency in a nuclear starburst, 
at the same time feeding the central BH at an Eddington-limited 
rate. This process consumes a large amount of the original cold 
gas reservoir. The combined energy generated from supernova 
explosions and the central engine wrecks havoc on the rest of 
the interstellar medium in the galaxy, shocking it to high tem- 
peratures and redistributing it to large scales, thereby shutting 
off further star formation and accretion. While the specific for- 
mulation of the problem and the methodology for solving it may 
vary from study to study (e.g., Granato et al. 2004; Springel et 
al. 2005; Croton et al. 2006; Hopkins et al. 2006; Hopkins & 
Hemquist 2006; Sijacki et al. 2007; Di Matteo et al. 2008), 
one generic prediction runs constant: the onset of nuclear ac- 
tivity, especially during the peak phase of accretion when the 
central object unveils itself as an optically visible AGN, in con- 
cert with supernova feedback from the accompanying starburst, 
liberates so much energy on such a short timescale that the bulk 
of the cold gas gets expelled from the galaxy. None of the ex- 
isting models makes very precise predictions about the cold gas 
content during the evolution of the system, but it seems clear, 
regardless of the details, that during the optically visible (unob- 
scured) phase of the AGN the host galaxy should be deficient in 
cold gas. This expectation is inescapable if AGN feedback is to 
have as dramatic an effect on the evolution of the host galaxy 
as has been suggested. 

Our observations present a challenge to the framework of 
AGN feedback just outlined. Far from being gas-deficient, the 
host galaxies of optically selected type 1 AGNs are, if any- 
thing, unusually gas-rich, and at the very least we can state 
with confidence that their gas content is normal. To be fair, 
the vast majority of our sample consists of Seyfert galaxies. 
They are hosted in disk (SO and spiral) galaxies, which proba- 
bly never experienced many major mergers, and their nuclei fall 
far short of luminosity threshold of quasars, which, although 
rare, dominate the BH mass density in the Universe (e.g., Yu 
& Tremaine 2002). These objects probably fall outside of the 



purview of major merger-driven models (Hopkins & Hernquist 
2006). Our sample, on the other hand, does contain seven ob- 
jects that formally satisfy the luminosity criterion of quasars, of 
which four (PG 0844H-349, PG 1426h-015, PG 2130h-099, and 
RX J0608.0H-3058) are detected in H I. All four are very gas- 
rich, with H I masses ranging from Mh i = 5 x lO'' to 2 x 10^" 
Mq and Mh i/Z^B.host values that are not obviously low. 

Beyond the present set of observations, two other fines of 
evidence contradict the notion that optically selected quasars 
are gas-deficient. Scoville et al. (2003) performed an unbi- 
ased search for CO emission from all z < 0.1 quasars from 
the ultraviolet-selected sample of Palomar-Green (PG) sources 
(Schmidt & Green 1983), finding that the majority of them 
(~75%) contain abundant molecular gas, ranging from ~ 
10'' to 10'° Mq; those that were undetected have upper lim- 
its consistent with the detections. Bertram et al. (2007) re- 
port very similar statistics from their recent CO study of low- 
luminosity quasars selected from the Hamburg/ESO survey. 
These molecular gas masses are comparable to, if not greater 
than, those typically found in mid- to late-type spirals. Ho's 
(2005a) compilation of CO observations of other PG quasars 
further reinforces this point, and even more extreme molecular 
gas masses (Mh, ~ lO'"- lO" Mq) have been detected in high- 
redshift quasars (Solomon & Vanden Bout 2005, and references 
therein). The Milky Way, for reference, has Mh, = 2 x lO' 
Mq (ScoviUe & Good 1989). Of the three PG quasars that 
overlap with our H I sample, for instance, PG 1426h-015 has 
Mh, = 6.5 X lO"* Mq, PG 2130H-099 has Mh^ = 3 x 10'' Mq, 
and PG 0844H-349 has an upper limit of Mh, = 1 x 10^ Mq. 
Combining this with the H I masses, all three objects have a 
total (H I -1-H2) neutral gas mass in excess of 10'° Mq. Yet 
another way to recognize that PG quasars have plenty of cold 
gas comes from noting that substantial amounts of dust have 
been detected in them. Haas et al. (2003) appraised the dust 
content in a sample of 5 1 PG quasars by combining far-infrared 
continuum observations along with millimeter and submillime- 
ter data. Their derived dust masses range from Mjust = 7 x 10^ 
to 3 X 10^ Mq, with a median value of 8.8 x 10^ Mq. For a 
Galactic molecular gas-to-dust ratio of 150, this corresponds 
to Mh2 = 1.3 X 10' Mq; for a higher gas-to-dust ratio of 600 
(Young & Scoville 1991), Mh, = 5.3 x lO'* Mq. These indi- 
rect estimates of molecular gas mass are consistent with those 
derived from CO observations. 

To graphically compare the gas masses of PG quasars side- 
by-side with those of the lower luminosity AGNs in our sample, 
we converted the molecular and dust masses into approximate 
H I masses. We adopt, for the sake of illustration, Mh2/Mh i = 
3, a value characteristic of early-type spirals (Young & Scoville 
1991), and Mu^/Mi^st = 600. We calculated BH masses for the 
sources using the H/3 line widths from Bororon & Green (1992) 
and optical continuum luminosities from the spectrophotome- 
try of Neugebauer et al. (1987), and then inferred fi-band host 
galaxy luminosities using the correlation between BH mass and 
bulge luminosity given in Kormendy & Gebhardt (2001), for 
simplicity assigning all the luminosity to the bulge. As sum- 
marized in the bottom panel of Figure 5, PG quasars span al- 
most 4 orders of magiutude in Mm /Lb, host, from objects as 
gas-poor as the most extreme elfipticals to as gas-rich as the 
most gas-dominated dwarf irregulars. Although this exercise 
becomes progressively more uncertain for high-redshift objects, 
we can attempt a similar order-of-magnitude calculation for 
the lune high-redshift quasars fisted in Ho (2007b), using BH 
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masses listed therein and the molecular gas masses given in 
Solomon & Vanden Bout (2005). Under the same assump- 
tions adopted for the PG quasars, we find an average value of 
log(MH i/LB.hosi) ~ -0.9, in the middle of the range observed 
for PG quasars (Fig. 5). In fact, this estimate, which assumes 
a local value for the ratio of BH mass to bulge luminosity, is 
most Ukely a lower limit because high-redshift quasars in gen- 
eral (Peng et al. 2006a, 2006b) and this subset of objects in 
particular (Shields et al. 2006; Ho 2007b) seem to exhibit a 
higher ratio (by a factor of a few) of BH mass to host galaxy 
mass. 

Given the lack of specific theoretical predictions, we can- 
not attempt a quantitative comparison of our results with AGN 
feedback models. Nonetheless, we wish to stress that the typ- 
ical unobscured quasar possesses a gas reservoir normal for its 
stellar content and does not appear to be depleted in cold gas. 
Whether this poses a difficulty or not for AGN feedback models 
deserves further consideration. 

5. FUTURE DIRECTIONS 

The pilot study presented here can be improved and extended 
in several ways. A number of the H I observations would ben- 
efit from longer integrations. Roughly 1/3 of our sample was 
undetected in the current experiment, and the upper Umits on 
their H I masses remain relatively high; nearly half of the non- 
detections have limits above Mu i = lO"' Mq. It would be use- 
ful to improve these limits by a factor of a few, down to Mu i 
Ri 5 X 10^ Mq, the scale of the Milky Way. Higher signal-to- 
noise ratio spectra of the detected sources would allow us to 
better study the H I line profiles, especially to quantify the ap- 
parent excess in line asymmetry and its implications for the dy- 
namical state of the gas and the circumgalactic environment of 
the host. In terms of new observations, it would be particularly 
worthwhile to observe a larger sample of high-Eddington ra- 
tio systems in order to verify the trends with Eddington ratio 
identified in this study. This can be easily accomplished us- 
ing the current SDSS database (Greene & Ho 2007b, 2007c). 
To put more stringent constraints on AGN feedback scenarios, 
it would be highly desirable to increase the number of high- 
luminosity AGNs (quasars) in the survey. Given the rarity of 
nearby quasars and the redshift limit of Arecibo, this will be a 
difficult task, but some progress can still be made. 

The current project represents a major initiative in character- 
izing the cold gas content in type 1 AGNs. We have demon- 
strated that H I observations can yield many useful insights into 
the nature of AGNs, BHs, and their host galaxies. This tool 
should be applied to other classes of AGNs. A useful extension 
of this program would be to observe a well-matched sample 



of type 2 (narrow-line) AGNs in order to test their relationship 
to type 1 objects. Kim et al. (2006; see also Ho 2005b and 
Lai & Ho 2007), for example, have argued that type 2 quasars 
are the evolutionary precursors, as opposed to simply misori- 
ented counterparts, of type 1 quasars. If true, type 2 sources 
are likely to be more gas-rich than their type 1 counterparts, or 
perhaps their interstellar medium might be less dynamically re- 
laxed, differences that can potentially be discerned through H I 
observations. 

The analysis presented here critically depends on a number 
of quantities derived from optical images of the host galaxy, 
namely its inclination angle, isophotal diameter, luminosity, 
morphological type, and, to a lesser extent, local environment. 
All of these measurements can be greatly improved with the 
help of better optical imaging. In many instances, the SDSS 
images (Ho et al. 2008) simply lack the necessary depth or 
resolution to yield clear-cut measurements of these parameters. 
High-resolution (e.g., Hubble Space Telescope) images would 
be particularly valuable, as they would simultaneously yield a 
reliable decomposition of the AGN core from the host galaxy, 
as well as a quantitative measurement of the bulge component. 
All of the correlations presented in this study should be reex- 
amined once such data become available, to see if their scatter 
can be reduced. 

We have stressed the utihty of H I observations, but, of 
course, a full inventory of the cold interstellar medium should 
include the molecular component as well, which, depending on 
the galaxy type, can dominate over the atomic component. A 
systematic CO survey of the objects in our H I program would 
be highly desirable, not only to complete the gas inventory but 
also to bootstrap the kinematical analysis to the CO line width. 
UnUke the current limitations of H I observations, the rotational 
transitions of CO can be observed to almost arbitrarily high red- 
shifts (Solomon & Vanden Bout 2005), a necessary ingredient 
to access AGN host galaxies at early epochs. Ho (2007b) dis- 
cusses the use of the CO Tully-Fisher relation as a powerful tool 
to investigate the host galaxies of high-redshift quasars. 
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